Geometric phases in semiconductor spin qubits: Manipulations and decoherence 
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We describe the efltect of geometric phases induced by either classical or quantum electric fields 
acting on single electron spins in quantum dots in the presence of spin-orbit coupling. On one 
hand, applied electric fields can be used to control the geometric phases, which allows performing 
quantum coherent spin manipulations without using high-frequency magnetic fields. On the other 
hand, fluctuating fields induce random geometric phases that lead to spin relaxation and dephasing, 
thus limiting the use of such spins as qubits. We estimate the decay rates due to piezoelectric 
phonons and conduction electrons in the circuit, both representing dominant electric noise sources 
with characteristically differing power spectra. 

PACS numbers: 72.25.Rb,71.70.Ej,03.65.Vf 



I. INTRODUCTION 

Recent demonstrations of coherent single-electron spin 
control and measurement in semiconductor quantum 
dotsii^i^ represent milestones on the way to quantum- 
state engineering with spin qubits^. In addition, work 
on coherent spin transport in nanostructures^'^'^ has 
revealed new possibilities for next-generation spintronic 
devices. The key behind these emerging technologies is 
the long spin coherence time in semiconductor materials. 

The standard technique for addressing and manipu- 
lating spins in semiconductors is electron spin resonance 
(ESR) controlled by external ac magnetic fields^. Alter- 
natively, effective internal magnetic fields can be gener- 
ated via the spin-orbit (SO) interaction by applied elec- 
tric fields. Proposals for coherent control of confined 
electron spins based on this combination have been put 
forward^iiSiiiiii^iH, and experimental progress has been 
reportedi^ii^. At the same time, SO interaction makes 
the electron spin sensitive to the electric noise ubiquitous 
in typical solid state environments^^iii. The combination 
of both provides an important mechanisms by which elec- 
tron spins decay and lose coherenc o-'^^i^^'^^'^^i^^'^'^ . Such 
electric field fluctuations are generated, e.g., by lattice 
vibrations, but in low magnetic fields the Nyquist noise 
in the electrodes may be dominanli^^. 

The details of the spin-orbit mediated interaction be- 
tween the electric field and the electron spin hide some 
interesting twists. As we discuss below, in suitable time- 
dependent electric fields the spin of an electron confined 
in a quantum dot acquires a non-Abelian geometric phase 
(a generalized Berry phase). On the one hand, this geo- 
metric quantum state evolution allows for new spin ma- 
nipulation strategies purely controlled by electric fields. 
In comparison to the alternative, ESR manipulation by 
ac magnetic fields, geometric spin manipulation is poten- 
tially more robust, since it is not affected by gate timing 
errors and certain control voltage inaccuracies. 

On the other hand, in fluctuating electric fields the 
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FIG. 1: Spin relaxation and decoherence rates in semicon- 
ducting quantum dots. Results of the present work (solid 
lines) are compared to those of Ref. [2l] (dashed lines), where 
the effect of piezoelectric phonons is analyzed to second order 
in the coupling to piezoelectric phonons. We find that geo- 
metric dephasing (evaluated to fourth order) leads to a satura- 
tion at low magnetic fields. Deviations between both results 
at intermediate field values are mainly due to the Nyquist 
noise of the conduction electrons. At higher magnetic fields 
(not shown) the results coincide. Parameters correspond to a 
parabolic quantum dot in GaAs with uiq = 5K at T = lOOmK. 



accumulation of random geometric phases leads to spin 
relaxation and decoherence for low and even vanish- 
ing magnetic fields^, predominantly induced by high- 
frequency noise, huj k^T ^ ^jlbB. This saturation of 
the spin decay rate at low fields has been overlooked in 
the literaturei^iSiiS^j although a similar connection had 
been discussed for free electrons in the presence of dis- 
order scattering^^. In comparison to the previously dis- 
cussed spin-orbit mediated mechanism o^^'^^ the geomet- 
ric spin decay is of higher order in the electric field. It re- 
quires a minimum of two independent noise sources cou- 
pled to two non-commuting components of the electron 
spin, whereby the non-Abelian character of spin rotations 
(properties of the SU(2) group) becomes relevant. 

The present paper is devoted to a detailed study of the 
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spin evolution of a confined electron in a time-dependent 
electromagnetic environment. It extends earlier work2^ 
which is based on a perturbative approach that allows 
handling quadrupolar and octupolar fluctuations but is 
restricted to small and adiabatic electric fields. Here 
we develop alternative adiabatic methods, that allow de- 
scribing large electric fields and corresponding displace- 
ments. 

To begin with, wc consider in Sec. Ill Al a "semi- 
classical" electron moving along a fixed trajectory. This 
limit describes the situation where an electron is confined 
in a narrow quantum dot potential which is shifted by ap- 
plied classical electric fields. The spin-orbit interaction, 
which is treated in leading order in the small parame- 
ter xo I Iso (dot size over SO length) , induces a geometric 
spin precession tied to the electron's motion. We illus- 
trate how electric fields can be used to manipulate the 
spin state via the SO interaction. We further provide 
a qualitative picture for geometric dephasing (Sec. Ill Bl 
and III C|) . if the fields are classical fluctuating fields com- 
pletely characterized by their power spectrum. A more 
rigorous analysis based on a fully quantum-mechanical 
treatment of the electromagnetic field follows in later sec- 
tions. 

Next, in Sec. IIII Al we treat the confined electron 
fully quantum-mechanically taking into account the spin- 
texture of the confined orbitals in a (parabolic) quantum 
dot. This approach allows us to study larger dots with 
intermediate and high values of xq/Iso- Spin-dressing ef- 
fects in larger dots lead to a renormalization of the g- 
factor that also influences the geometric spin evolution. 
We first consider the effect of a classical applied field and 
then generalize to the case where also the electromagnetic 
field is treated quantum-mechanically within a path in- 
tegral approach. This method, presented in Sec. IIIIBI 
allows us to compute the dephasing and relaxation rates 
of the electron spin by means of a systematic diagram- 
matic perturbation theory. 

In section ITVl we present our results for the spin relax- 
ation and decoherence rates in semiconducting quantum 
dots. We find that at low temperatures and low magnetic 
fields, B < IT, Ohmic fluctuations originating from the 
electrodes dominate over the phonon effects in their in- 
fuence on the spin relaxation^. This is simply due to 
the fact that the spectrum of excitations of a metallic 
electrode is much denser than that of phonons at low en- 
ergies. At still lower magnetic fields and temperatures, 
however, the dominant mechanism leading to relaxation 
is the accumulation of random geometric phases. This 
leads to a saturation of the decay rates e.t B — 0. Both 
features are included in the data shown in Fig.[T]and com- 
pared to previous predictions. Further results including 
the temperature dependence under typical experimental 
conditions are presented in section HVl 

Finally, in Sec. |Vl we study an electron in an array of 
quantum dots and show that in the presence of SO inter- 
action arbitrary spin rotations can be reached by a series 
of coherent electron tunneling processes in a multi-dot 



geometry. We discuss the possibility to use such multi- 
dot systems for geometric electron spin manipulations. 

II. GEOMETRICAL SPIN EVOLUTION IN A 
SEMICLASSICAL PICTURE 

A. Geometric phases and spin rotations 

Geometrical phases, introduced in the pioneering 
work of Berry'^?., describe how a quantum system in a 
non-degenerate state evolves when adiabatically driven 
around a closed loop in a control parameter space. Apart 
from a dynamical phase, depending on the energy and 
elapsed time, the state acquires a contribution, called 
the Berry phase, that only depends on the geometry of 
the loop in parameter space. An important extension of 
the Berry phase covers the case when the initial state 
belongs to a degenerate subspace^^. Then the geometric 
phase is replaced by a non-Abelian unitary transforma- 
tion of the initial state within the degenerate subspace, 
and again, this unitary transformation depends only on 
the geometry of the loop. 

This non-Abelian generalization of the Berry phase is 
of direct relevance for an electron confined in a quantum 
dot: Time reversal symmetry in the absence of a mag- 
netic field implies that the spectrum of a quantum dot 
is a collection of time-reversed Kramers doublets even in 
the presence of spin-orbit interaction. In the absence of 
SO coupling these states correspond to pure spin states. 
In the general case the doublets persist, and we call them 
pseudospin states. 

The adiabatic variation of the position (and possibly 
shape) of the potential confining the electron is described 
by a time-dependent Hamiltonian H{t). Adiabaticity 
ensures that an electron which at time i = is in the 
ground state doublet remains within the corresponding 
subspace of instantaneous ground states of H[t). The 
adiabatic time evolution thus appears as a non-Abelian 
SU(2) transformation within the ground state doublet^. 
In the presence of spin-orbit coupling the unitary trans- 
formation results in a geometric spin rotation, which only 
depends on the trajectory of the electron in real space. 

To illustrate the geometric spin rotation in more detail 
we consider a 2DEG, assumed to be grown along the [001] 
direction. An electron in the 2DEG experiences a spin- 
orbit coupling which takes the form {h = 1, x axis along 
the [100] direction) 

Hso = a [Pydx - Pxd-y) + (3 {PyCfy " Px^x) = —pX^o^ ■ 

(1) 

Here p denotes the momentum of the electron in the xy 
plane of the 2DEG and (t/2 is the spin operator. The 
Rashba (a) and linear Dresselhaus (/3) couplings can be 
lumped into the spin-orbit tensor Ago 

A.-'.,„(-/7), (2) 



An implicit summation over the x, y coordinates is as- 
sumed, i.e. pK^^o- = J2^,u=x,yPt^ (Ko^),,^^'^- The ten- 
so r Ago set s the s cale for t he spin-orbit length 4o = 
a/ I det Asol = {m^y\a'^ — that defines the typical 

length scale of spin textures. In typical GaAs/AlGaAs 
semiconductor heterostructures we have 4o ~ 1 ~ bfiuiM 

We further assume that by lateral structuring of the 
2DEG a parabolic quantum dot, V{f) — ^muj^r'^ , is 
created with energy scale ujq related to the orbital size 
~ Xj ^mujQ. Typical dot sizes are in the range 
~ 30 — lOOnm. Hence xo/lso is usually small, of the 
order of 0.1 — 0.01. The electron is assumed to be strongly 
confined in the ground state of the potential. By apply- 
ing electric fields (modifying the confining potential) we 
move the dot and the electron along a trajectory Rc{t). 

In the considered limit, xq/Iso ^ 1, the degenerate 
doublet of states of the confined electron simply dif- 
fer in the spin orientation. Eq. ([1]) then implies that 
the confined spin experiences an effective magnetic field, 
Hso = Bso{t) ■ 6-, with Bso = -p^{p)\~a^. For a strongly 
confining potential we have (p) ~ mRc. Therefore the 
induced spin-evolution is described by"^^ 



C/ad = T exp i-i dt Rc^so ^ 



P exp ^— i y dRcXso ^ 



(3) 



with T and P denoting time- and path-ordering opera- 
tors, respectively. The tag 'ad' refers to the constraint 
of adiabatically slow paths, \tic\ ^ a^owo, which guar- 
antees that the evolution affects the electron spin but 
does not induce transitions to excited dot states. Clearly, 
?7ad depends only on the geometry of a given path C it- 
self, not on the time dependence of Rq- As we shall 
show in Sec. IIII A| this result can be generalized to larger 
parabolic dots for which spin texture effects are impor- 
tant. The main difference in that case is that A^^^ is 

replaced by a renormalized tensor A^^ . 

In order to visualize the spin rotation described by 
Eq. ([31) it is instructive to use the connection between 
SU(2) and S0(3) rotations. The change of the orienta- 
tion of a sphere rolling on a plane along a path Re is 
characterized by a rotation operator 
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(4) 
(5) 



where the components of A denote the standard SO (3) 
generators. Comparing ([3]) and ^ and choosing the 
radius of the sphere to be Rq = Zso/2 we note that 
a spin rotation C7ad associated with a path Rq maps 
onto a rotation of the sphere rolling along a trajectory 
Rc'{t) = 2Rc{t)Xsph (see Fig. H]). If only Rashba cou- 
pling is present then the paths C and C have the same 
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FIG. 2: (a) The geometric spin precession due to SO interac- 
tion for an electron adiabatically moving within a 2DEG (left) 
is equivalent to the changing orientation of a sphere rolling on 
a plane (right) . Both paths are related in a one-to-one fashion 
to each other, (b) The length of a straight path required to 
perform a spin flip versus the angle of the path relative to a 
crystal axis for several ratios (3/ a with a constant Iso- The 
area of the ellipses remains constant and equal to irlso- 



shape. In general, if also Dresselhaus coupling is present, 
then C and C have different shapes, but the qualitative 
analogy persists. 

In a 2DEG, with both Rashba and Dresselhaus cou- 
plings present, the change of the spin orientation induced 
by moving the confined electron along a straight line de- 
pends on the direction of trajectory relative to a crys- 
tal axis. The distance to perform a spin-flip is strongly 
anisotropic for a ^ f3, being enhanced in the [110] direc- 
tion. This property is illustrated in Fig.[2l^b). 



B. Spin manipulation and spin relaixation 

To illustrate how the geometric phase evolution given 
by Eq. Q can be used to manipulate the spin state we 
consider a quantum dot orbiting N times around a closed 
circular path, R(t) — Rq{cos cut, sin ut^O), with period 
to — 2tt/ll}. The spin-orbit length is assumed much larger 
than the loop size, Igo 3> Rq- For simplicity we assume 
pure Dresselhaus coupling. The resulting spin precession 
is then described by Uad = [Uad{to)]^ , where the evolu- 
tion operator of a complete loop is 



C/ad(to) = P exp 



27r 



Rq 



i<7x sin ( 



CT„ cos ( 



To second order in Rq/Iso it reduces to 



C/ad(to) 



1 — ia 



2ttRI 



exp 



'2 



(6) 
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where A{to) = ttRq denotes the area enclosed by a loop. 
As will be shown in Sec. IIII Ai Eq. ^ is valid for small 
closed loops of arbitrary shape. The full time evolution 
after N circles can be approximated as 

i 4:A{t) . 



ad 



■ exp 



2 12„ 



(7) 



where A{t) — N A(to) is the area swept hy N — t / rep- 
etitions of the loop. Obviously the motion in small closed 
loops corresponds to an effective static magnetic field, 
Bcs = 4A(to)/(io ^so) pointing along the z-direction. 

Eq. ^ also holds if the trajectory Rc{t) is a stochastic 
variable, typically induced by a fluctuating electric field, 
Rc{t) cx E. In this case the area A{t) in Eq. ([7]) is a 
random variable whose dispersion increases linearly in 
time, (A(t)2) « Aq t. Assuming that A{t) is Gaussian 
distributed, we conclude that off-diagonal elements of the 
spin density matrix decay as 

Am' 



exp 



(8) 



Thus the accumulation of - random - geometric phases 
leads to dephasing, even in the absence of external mag- 
netic fields. The coefficient Aq introduced above depends 
on the amplitude and frequency of typical electromag- 
netic fluctuations, and its proper computation requires a 
fully quantum-mechanical treatment of the electromag- 
netic field^, reviewed in Section fill Bl At this point we 
note that the area itself is proportional to the square of 
the electric field, E'^ , with oj ^ T the typical frequency of 
electromagnetic fiuctuations, and the characteristic time 
for making a closed loop is ~ l/w. Hence, the geomet- 
rical relaxation rate is proportional to (w E.g., 
for Ohmic fluctuations we have \E^\'^ ~ and conse- 
quently the corresponding relaxation rate scales as T^. 
Note that this geometric dephasing is of fourth order in E 
and, therefore, requires fourth order perturbation theory 
in the electromagnetic field. 



C. Spin relaxation in a magnetic field 

So far we assumed that no external magnetic field is 
applied. Now we study the relaxation and dephasing 
of a spin in an in-plane magnetic field B in the pres- 
ence of a classical stochastic field, Bso — RcK^o — 
(Bsox, Bsoy, Bsoz), characterized by the power spectrum 



Sij{uj) 



dte^"*(B,o,(i)5soj(0)). 



(9) 



The external field B induces a precession of the spin 
around its direction with Larmor frequency ujb = 
—giXB\B\, while the fluctuating fleld B^o leads to relax- 
ation and dephasing with rates Fi = T^^ and r2 = Tj"^, 
given to leading order by the standard expressions22iSi^ 



(2) _ 



.(2) 



(2) 



2^ll(0) . 



(10) 
(11) 



Here 5'_L(a;) and 5*11(0;) are the spectral densities of the 
components of B^o perpendicular and parallel to S, 
and the superscript (2), introduced to distinguish from 
later extensions, refers to second order. In case of 
isotropic random motion of the dot and purely Rashba 
(or purely Dresselhaus) spin-orbit coupling we have 
S^iuj) = 5||(w) = S{uj). Typically S{u} = 0) = if the 
electron remains localized in space, since only unbounded 
paths can produce a static Bso Rc ■ Therefore the last 
term of Eq. pip vanishes in this case and we obtain in 
leading order^ T2 = 2Ti. 

We can also estimate the relaxation rate induced by the 
Berry phase mechanism, which is of higher order than 
expressions Eqs. ([TU]) and pT|) : The effective magnetic 
field generated by the Berry phase term is proportional 
to i?off ^ A ^ Rc X Rq, and is perpendicular to the 
in-plane magnetic field. The relaxation rate of the spin 
to fourth order in the coupling to the stochastic field is 
proportional to the Fourier component of the autocorre- 
lation function of this effective field at frequency w ~ ws. 



(12) 

where C\\ and refer to the spectral functions of the 
components Rc,\\ £ind Rc,±, parallel and perpendicular 
to the in-plane magnetic field. 

Strictly speaking, Eqs. (|10m2p hold only in the limit 
lub ^ Fi,r2. In the opposite case, ujb < Fi,F2, also 
referred to as the Zeno regime, the relaxation times Ti 
and T2 lose their meaning, and other dissipative rates 
need to be considered. Although r\ and F2 vanish fast 

in the limit ujb 0, F^^-* saturates and scales to a finite 
value. This is because the Berry phase relaxation can also 
be induced by independent fluctuations of frequency uj ^ 
T ^ UJB, that have a 'beating' (frequency mismatch) at 
frequency ws, which is in resonance with the Larmor spin 
precession. As a result, for sufficiently small ujb but finite 
T one always ends up in the Zeno regime, dominated by 
the Berry phase term, Eq. ([l3 



III. 



BEYOND THE SEMICLASSICAL 
APPROACH 



In this section we provide a fully quantum-mechanical 
treatment of a confined electron's spin. We do this in two 
steps: First, we describe quantum-mechanically the spin 
and orbital state of the confined electron in a classical 
time-dependent electric field. Then, we generalize this 
within a path integral formalism to include quantum- 
mechanical fluctuations of the electric field. 
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A. Adiabatic approach in a classical electric field 

1. Motion in a parabolic dot 

To be specific, let us first consider an electron in 
a parabolic confining potential V{r) = ^rruJ^r'^ with 
energy scale loq corresponding to a typical orbital size 
= 1/ ^muj^. For the moment we will assume that 
the magnetic field is zero. In the presence of an ap- 
plied or fluctuating classical in-plane electric field E{t) 
the Hamiltonian takes the form 

n{t) = ^ + V{f) + eE{t)r + nso 
Am 

= ^ + V{r-Rc{t))+H,o+C{t) , (13) 

with C{t) ~ E'^it) being a time-dependent constant of no 
relevance, and Rc{t) = —eE{t)/mujQ. Thus the effect of 
the homogeneous electric field is to move the center of 
the potential along a trajectory, Rc{t). We assume that 
-Rc(O) = 0. Then the Hamiltonians for times t > are 
related by the displacement operator W{t) = e^vRcit)^ 

H{t) = Wit)HiO)W+{t) + Cit). (14) 

The center of the potential, Rc{t), is our external con- 
trol parameter in the language of geometric phases. 
Eq. allows us to construct the instantaneous eigen- 
states {|T„(t))} of Hit) from the eigenstates {|t„)} for 
Rc = 0, as {\Tn{t))} = {>V(<)|t„)}. The states come in 
Kramers doublets with n = 0, 1, 2, . . . and r =| , | denot- 
ing their pseudospin orientation. We now introduce the 
'adiabatic states' or 'center of motion states' as 

\^it))^w+im{t)) . (15) 

Here the state |4'(i)) satisfies the time-dependent 
Schrodinger equation with the Hamiltonian Ti.{t), while 
the adiabatic states |^'(i)) move together with the dot 
and satisfy the equation of motion 

idt\^{t)) = {m^+{t)w{t)+n{o)+c{t))\^{t)) . (16) 

Using the explicit form of the operator W we thus find, 
apart from a trivial overall phase generated by C{t), that 
the evolution of \^{t)) is described by the effective Hamil- 
tonian 

Hee=n{0)-Rcit)p, (17) 
and the corresponding evolution operator, 

U{t)^Tcxpi^-zJ^ dt'[n{0)-^-p]^ . (18) 

Clearly, from the above equations it follows that if at 
time t — the electron occupies the ground state 
Kramers doublet, i.e., |*(t = 0)) = Yl^T-a^it = 0)|to), 
and the external perturbation changes sufficiently slowly 



in time, |J2c(i)| <C XqUq, then the second term in 
Eq. p?|) cannot generate transitions to the excited states 
of 7^(0), and |^(i)) stays within the ground state dou- 
blet, \^{t)) « X^r '^i"(*)l'^o)- Under these conditions the 
effective Hamiltonian can be approximated as 

H,s~~Rc{t)PopPo, (19) 

where Pq = J2t I'^o)(7'o| is the projector to the ground 
doublet of H{0), and we took the energy of the ground 
state to be Eq = 0. The corresponding adiabatic evolu- 
tion operator within the center of motion ground state 
subspace then reads 

U^it) = Texp (^ij^ dt'Rc{t)PopPo^ ■ (20) 

Although not obvious, the spin-orbit coupling plays 
an essential role in Eqs. and ((^D|). As shown in 

Appendix \X\ we have, 

PopPo = -X-^PoaPo - A-i^f, (21) 

where we introduced the non-trivial 2x2 spin-dressing 
tensor-- z, that relates the matrix elements of the spin 
operator to the pseudospin operator t. Thus an adia- 
batic motion induces a pseudospin rotation within the 
ground state multiplet. Note that the formulas above 
are not perturbative in the spin-orbit coupling, and they 
hold for parabolic dots of any size compared to the spin- 
orbit length as long as the external field fiuctuations are 
adiabatic. Thus the effects of strong spin-orbit coupling 
can be fully taken into account by replacing by the 

dressed spin-orbit tensor A^^ = A^^z, so that the adia- 
batic evolution operator becomes 

i/ad = P exp £ dRc\,'f^ , (22) 

with P being again the path ordering operator. From 
Eq. it is obvious that the renormalization of the 

g-factor is also reflected in the renormalization of the rel- 
evant spin-orbit length scales, Zso — > 4o = \/ \ det Xso \ ■ 

This renormalized Zgo has been shown in Fig. [3] as a 
function of the ratio xq/Iso- In small quantum dots, 

xq/Iso ^ 1, one has ~ 5^^, and therefore A^^ = A~^, 
in agreement with the assumptions of Sec. [Ill 

2. Connection with perturbation theory 

The preceding discussion relied on the specific model of 
a parabolic quantum dot, which enabled us to construct 
the instantaneous eigenstates of the Hamiltonian for arbi- 
trarily large but slowly changing displacements, Rc ~ E. 
An alternative approach to construct C/ad that is valid 
for an arbitrary shape of the confining potential is to use 
perturbation theory in the driving term iJjnt = eE{t)r 
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FIG. 3: (Color online) The dressed spin-orbit length Iso as a 
function of the typical orbital size xo — \J 1 /muo for different 
ratios /3/a. Exchanging a and /3 yields identical curves 



of Eq. to determine the instantaneous eigenstates of 
H{t)^^. This perturbative approach has the disadvantage 
of breaking down for displacements Rc comparable to the 
typical dot size Xq. However, in a conventional setup of 
single quantum dots defined by lithographic gates, one 
is usually constrained to small displacements \Rc\ 
such that the perturbative approach remains useful. 

In the absence of external magnetic fields the pertur- 
bation iJint = eE{t)r does not connect states within the 
same doublet due to time reversal symmetry. One can 
then approximate the lowest energy instantaneous eigen- 
state \To{t)) as 



\Mt)) 

Wpcrt(i) 



J2{\ro){ro\+eE,it)J2 



where for brevity we have denoted by En the energy 
of doublet n at time t = 0, and |t„) = |t„(0)). Then 
the effective adiabatic Hamiltonian, describing the evo- 
lution of the center of motion wave function |4'(t))pert = 
Wport(^)l*(^)) within the ground state multiplet, reads 



n>0 ^ 



xPnV - yPnX 



Eo) 



Pn 



where the operators P„ project to the n'th doublet state 
of the unperturbed Hamiltonian 7i(0). Furthermore, in 
this equation we dropped full differential terms involving 
Efj,Ej, + Ef^Ejj that give no contribution for closed paths, 
and we denoted hy x,y the components of the electron 
position operator. The above expression coincides with 
the B — Q limit of the results obtained previously using a 
more general perturbative approach in Ref. '2^.— Treat- 
ing the general case, B ^ Q, within this approach requires 
some care since the appropriate perturbation theory is 
quasi-degenerate22i2^. 

While it is not obvious, the above perturbation expan- 
sion approach and the displacement operator results pre- 
sented in the previous subsection are equivalent. Due to 



the slightly different choice of the transformation yV+ [t) 
the two approaches yield different adiabatic evolution op- 
erators, however, both Hamiltonians describe the same 
physical pseudospin precession, and physical observables 
are independent of this choice of basis. To illustrate this, 
let us compare the two matrices U^d after a path that 
starts and ends at the origin. For a closed path Rc{t), 
the perturbative approach yields a pseudospin transfor- 
mation that can be expressed as 



P exp 



2 Jc 



{dRc X Rc). CP"* 



^pcrt . ^ , _ „-^^ \^ 



n>0 



{TQ\xPnV - y^n^ko) 

{En - Eo)-' 



(24) 



As we show in Appendix |^ for the particular case of 
a parabolic dot sum rules imply that C^"' simplifies to 

C^y''* • o-rr' = -ixlirolp^PoPy - pyPaPxlrl^) (25) 

On the other hand, the exact displacement operator 
approach of the previous subsection yields 



P exp 



dRc ■ PopPo 



(26) 



For a closed path with a small enclosed area we can ap- 
proximate this expression by expanding to second order 
in the exponent and re-exponentiating as 



U^d = P exp \i [{dRc X R, 



ic)z X 



(27) 



X 2 {PoPxPoPyPo - PoPyPaPxPo) \ ■ (28) 

Eqs. ISll), ([281) and Eq. ([251) imply that both approaches 
yield identical results. Furthermore, in the case of a 
parabolic dot we can show using the alternative form for 
[/ad, Eq. (P^ . that for small closed paths the evolution 
operator J7ad of the unshifted states, |^(i)) is approxi- 
mately given as 



ad 



ad 



exp 



i AA 
2F"^ 



(29) 



where A — ^ ^^^{dRc x Rc)z is displacement area 
spanned by the path. This generalizes the semiclassi- 
cal result of the previous section to the case of arbitrary 
spin-obit coupling. 



B. Treating the electromagnetic field 
quantum- mechanically 

In the first part of this section we assumed that the 
time-dependent electric field acting on the electron is 
classical. However, if the typical frequency of the field 
reaches the temperature then the electric field must be 
treated quantum-mechanically. To treat the quantum 
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fluctuations of the electromagnetic field, we shall employ 
a hybrid formalism, where we describe the fluctuations 
of the electromagnetic field within the path integral ap- 
proach, while the quantum dot shall be treated within 
an operator formalism. Here we only summarize the re- 
sults, the details of the calculation are presented in Ap- 
pendix [b] 

For simplicity, we consider an electron spin in a 
parabolic quantum dot and dipolar electric fields as in 
subsection IIII A 11 Then the adiabatic evolution opera- 
tor that describes the evolution of the dot state for 
given initial Ei and final E / states of the electromagnetic 
environment within the ground doublet is given by 



Ef\U,d{t)\E, 



(30) 



w 



Hz - -^(IbBPo&Po 



(31) 



Hait) = Rc A-iPoo-Po = -^A-^Pq^Po • (32) 

Here the functional integral is performed over all pos- 
sible fluctuations of the bath, E{t), each of them cor- 
responding to a different path Rc{t) between times 
and t. We also allowed for the presence of an in-plane 
magnetic field S, small compared to the level spacing 
of the dot, ~ Wo- The adiabatic displacement opera- 



tors W = e'P^-^/'^ 



_ p~ipRc 



in Eq. (|5D|) transform 



the wave function of the dot to the center of motion ba- 
sis, l^'(i)) — > \^{t)), where the evolution is described by 
Eq. p^ . projected by the projector Pq to the lowest lying 
two eigenstates of 7i(0), in the absence of the magnetic 
field. 

Thus the pseudospin evolution under a quantum elec- 
tric field is the coherent sum of classical evolutions over 
all possible field fluctuations weighted by the action of 
the decoupled bath, S'^. We remark here that S'g is not 
the non-interacting bath action, Sb, but contains a cor- 
rection AHb = —e^\E{t)\'^ / (2to(jJq) due to the back reac- 
tion of the dot (see Appendix [B|) . For a Gaussian S'^ we 
thus mapped the evolution of the spin to the well-studied 
spin-boson model, where, however, external bosonic fluc- 
tuations are coupled both to the x and y components of 
the spini^^i^ 

To describe dephasing and spin relaxation at the fully 
quantum-mechanical level we generalized the formulas 
above to the evolution of the reduced density matrix for 
the pseudospin in the center of motion basis, 



PD{t)-i 



(rolTri 



yv+p(t)yv 



(33) 



where Tr^ stands for a trace over the bath degrees of 
freedom, p{t) = U{t)p{0)U^{t) is the density matrix of 
the complete dot-bath system, and W = exp(— ip • Rc)- 
Note that now Rc — ~eE/mujQ is an operator instead of 
a c-number, and that W = exp(— ip • Rc) acts both on 



the bath and on the quantum dot. Then spin relaxation 
and decay corresponds to the diagonal and off-diagonal 
components of pD{t)TT', respectively. The initial state 
p{0) of the system should not affect the dynamics at 
long times'^'* , so we will assume the dot-plus-bath system 
to start at t = in a well defined ground doublet and 
the bath in thermal equilibrium, p{0) = >V^/o(0)yV = 
Pd{0) <E) pb{0)- This choice is technically convenient and 
physically describes a definite initial state with slightly 
entangled bath and dot states. Then the evolution of 
PD{t) can be described as 



Mt) 



T 



K 



Pd{0) 



(34) 



where now z runs along the usual double branch time 
contour K running from — to t and then back to 0, 
as depicted in Fig. [7] in Appendix [Cl and Tr- is time 
ordering along this contour. Here (. denotes the 
average over the electromagnetic field fluctuations along 
the Keldysh contour, 



V[E] (E+{0)\pBmE-m e 



where E+ and E^ denote the electric field along the 
upper and lower branches of the contour and satisfy the 
boundary condition, E+{t) = E-{t). 

Unfortunately, one cannot integrate out the electric 
field exactly. However, one can use the systematic dia- 
grammatic approach of Ref. [1^ to do perturbation the- 
ory in Hcit) and obtain the relaxation rates within a 
Markovian approximation."*" The results of this calcula- 
tion are presented in the following Section. 



IV. RESULTS FOR RELAXATION AND 
DEPHASING 



We are now ready to summarize the results that are 
obtained by the formalism developed in the previous sec- 
tion. Some technical details of the calculation are pre- 
sented in Appendix [Cl We remark here the present sys- 
tematic approach confirms what had been derived within 
the perturbative formalism developed in Ref. 23. Unlike 
in said work, where numerical results for Ti were com- 
puted in a limited set of cases, the results for Ti and T2 
presented in this section are fully analytical, which there- 
fore allows for an explicit analysis of the decay rate de- 
pendence with all relevant physical parameters, such as 
magnetic field orientation, spin-orbit coupling or quan- 
tum dot size. 

We have computed the spin relaxation (Tj~*) and de- 
phasing {T2^) rates of confined electrons subject to an 
in-plane (renormalized) magnetic field B at an angle 9 
with respect to direction [100]. To fourth order in the 
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spin-orbit coupling we obtain 

^ = 2(ma;o)2(a2+/32-2a/3sin20Kcoth-^ 



A{ujb) 



+ 2{mxof {a'^ + (3Y + 4:a^P'^ cos 46 F+{uob) 

+ 2(mxo)*(a' - l3''fF-{ujB) , (35) 



To 



2Ti 



+ {mxaTia^ + /3 - 2a/?sin26')^i^(wB) . (36) 



Here ub — xq is the quantum dot size (as- 

sumed parabolic) and the functions and F are defined 



F+{iOB) 
F-{ujb) 
F{cjb) 



duj A I 



Ai 



-1 _ cosh{Q/2KBT) 
cosh{uJB/2KBT} 



dci,^, A{^ 

Stt'^ I _ cosh{0/2KBT) 



xRe 



— (I)''A(6:))^csch^ 



1 



Stt 



cosh(wB/2Ki3T) 

{2Kbt) 



1 



— — iO+ lub + iu + iO+ 



with the spectral function A{uj) of the dimensionless 
electric field, exoE/uo, defined in Appendix [Cl Func- 
tions F^ relate to the F^2 used in Appendix O by 
F± = {F^±F^)/2. Note that both F+{ujb) and ^(wij) 
vanish for tj^ 0, while -F~(0) remains finite. There- 
fore, this latter part of the fourth order contribution can 
be identified as the Berry phase contribution. Indeed, 
this expression is identical to the one we obtained previ- 
ously in Ref. [23|. Note also that this term is isotropic 
and does not depend on the direction of the magnetic 
field, while all other second and fourth order contribu- 
tions do so for a, 7^ 0. In the formulas above we 
assumed that the spectral functions of the components 
of the dimensionless electric field ex^E /ljq are isotropic, 
A^{lo) = Ay{Lo) = A{uj). 

The spectral function A{u;) gives the density of electro- 
magnetic excitations that contribute to dephasing, and it 
depends on the specific source of electromagnetic fluctua- 
tions. Important sources of electric field fluctuations are 
piezoelectric phonons^^, and at low magnetic fields Ohmic 
charge fluctuations^^, 

A(w) = Aph(w) + Ao(w). 

For the case of piezoelectric phonons an estimate 
of the spectral function Api-^{uj) of the induced dimen- 
sionless electric field is outlined in Appendix |DJ For 
GaAs/GaAlAs heterostructures we obtain Aph(w) = 
w^Aph^o/wo, with Aph — 2.5 • 10~^K^^nm^^. At higher 
magnetic fields (frequencies) these fluctuations give the 
dominant contribution to spin relaxation/dephasing. 

At low magnetic fields (low frequencies), on the other 
hand, Ohmic charge fiuctuations of the electric environ- 
ment of the 2DEG near the dot are expected to domi- 
nate. For these, the spectral function is Aq{uj) = Xqoj/ujq 



for Fermi liquid leads, with Aq ~ (e^//i)Re[Z], Z be- 
ing the impedance of the leads. It is rather difficult to 
compute the exact value of An, since it depends on the 
precise geometry of the leads and one must also take 
into account how the equilibrium electric field fluctua- 
tions inside the 2DEG extend to the quantum dot, and 
to what extent those fluctuations can be screened. It is 
not unreasonable, however, for typical sheet resistances 
of 10^ - lO^n/D to assume \n 10"^ - 10"''. 

In our calculations we have neglected the feedback cor- 
rection — e^|£^p/(2mwQ) that changes the action Sb 
S'b- In the case of the phonon bath, the effect of this 
term is to make phonons somewhat softer close to the 
dot. However, this polaron-type effect should be small, 
since the mass of the electron is negligible compared to 
the atomic masses. The feedback correction might be 
more important for the Ohmic bath, but will only affect 
An, not the Ohmic character of the bath. 



Implicit in the derivation of Eqs. (|35|) and (|36|) is the 
assumption that the effect of the thermal bath leads to 
an adiabatic evolution of the spin in the quantum dot. In 
physical terms one can anticipate that this implies a con- 
dition for the bath temperature fcsT <C wo, so that the 
dot is not heated above the ground doublet by the bath 
fluctuations. More precisely, the classical adiabatic con- 
dition Rc/xo ^ ujQ given in Sec. Ill Al translates in this 
context into the relation {kBT)^A{kBT) <C uJq. For typ- 
ical parameters in the case of the ohmic bath this means 
indeed 0.2kBT <C ujq. For the piezoelectric phonon bath 
this implies 6 (T[K])^/^ fcisT < luq, where T[K] is the 
bath temperature in Kelvin. Non-adiabatic corrections 
are beyond the scope of this work, although they could 
in principle be taken into account in the calculation of 
Appendix [Cl 



The flrst (second order) term in Eq. ((551) is propor- 
tional to ujbA{u!b) at small Larmor frequencies. Con- 
centrating on this term one concludes that the relaxation 
and dephasing rates vanish for cub q18j24j25_ ^pj-^g 
ond and third terms in Eq. (|36p . which are of fourth 
order in the spin-orbit couplings, involve a convolution 
of the spectral function at frequencies up to the temper- 
ature scale, and F^{lob) does not vanish at low fields and 
non-zero temperature. Therefore the relaxation rates sat- 
urate both for the phonon and the Ohmic baths as the 
field is lowered. We remark that the sum of the terms 
F+ and F^ resemble the formula, Eq. p^ . anticipated 
in Sec. Ill C[ although the full quantum treatment given 
here is necessary to identify precisely the functions C±_ 
and C|| in Eq. ()12l) and to arrive at consistent quantita- 
tive predictions. 

In Fig. m we plot the relaxation rates induced by the 
Ohmic and phonon baths separately for two different sets 
of parameters. For small magnetic fields huB kBT 
(but still above the Zeno regime) the relaxation rates 
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FIG. 4: Relaxation rates due to piezoelectric phonons and 
Ohmic fluctuations in two scenarios, both as a function of the 
magnetic field B applied along the [100] direction. For the 
strengths of Rashba and Dresselhaus spin-orbit coupling a ra- 
tio q://3 = 4 is assumed^. In (a) we plot the unfavorable sce- 
nario for quantum information precessing, with An = 5 - 10~^, 
Iso ~ 1500nm and a large dot xq ~ 115nm {ujo ~ IK). In 
(b) we plot the rates for favorable conditions An — 10"'*, 
Iso = SOOOnm and xq = 36nm (tJo = lOK). Lines range 
from a bath temperature of lOOmK (lower rates) to 900mK 
(higher rates) in steps of 200mK. Shaded in gray is the region 
Tj~^ > LOB ~ —gfJ.B\B\ for which the rotating wave approxi- 
mation fails and the Zeno regime sets in (see Sec. Ill Cp . Note 
the extreme dependence of relaxation rates on the specific 
conditions. 



"cr 10° 




« 10' 



10' 10' 10° 
COj/cOo 

FIG. 5: Total relaxation and decoherence rates for a 51nm dot 
{ujo = 5K) at T = lOOmK. Other parameters are Lo = 3/im, 
a/P = 4, e» = 0, An = 10"^ Three regimes are clearly visible. 
In the inset we give a polar plot of the dependence of the 
rates with magnetic field angle 9. This angular dependence 
is negligible in the saturated regime, and of the form cx 
a'^ + f3'^ - 2af3 sin 261 at higher fields. 




10' 10' 



FIG. 6: (Color online) Phase diagram of the three relaxation 
regimes at An = 10~^, a//3 = 4, /so = 3/im and ^ = 0. A win- 
dow of ohmic fluctuations - dominated relaxation and dephas- 
ing opens up at lower temperatures. The saturation regime 
in this plot is dominated by phonon fluctuations, although 
regions exist (not shown) at lower fields and temperatures 
where the saturation is mainly due to Ohmic fluctuations. 



reach a clear saturation regime, with values 



form characteristic of an Ohmic bath. 



{ksTf 



(37) 



numerical 
16.5, Cg 
2.58- 10". 



constants equal equal C4 
= 327rV21 ~ 466 and Cg 



where the 
87r3/15 ~ 
1287rVl5 ~ 

As the field is increased, the second order terms begin 
to dominate. At low enough temperatures Ohmic fluc- 
tuations become dominant first, since their weight grows 
faster at low fields than the phonon contribution. Thus 
below a certain temperature there is a window of mag- 
netic fields in which the relaxation rates have a scaling 



Tf 1 w 2T.2'^ cx uj% coth 



2kBT 



At still larger values of the magnetic field the phonon 
fluctuations dominate the second order relaxation chan- 
nel, and the rates assume a 'phononic' scaling form. 
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2kBT 



These cross-overs are shown in Fig. [5] for both relaxation 
and dephasing rates. The different regimes are summa- 
rized in a 'phase diagram', which we show for typical 
parameters in Fig. [Sj 

As noticed in Ref. [U, T^^ = 2X2^ up to second order 
in the coupling xq/Iso (i-e. at high fields). This relation 
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is violated by the fourth order terms due to the F con- 
tribution in (|36p . but is again restored in the saturation 
regime, where the fourth order term ^ F not considered 
in Ref. [2l| vanishes. 

For a, P second order terms that dominate the 
high field behavior have a strong dependence on 9: T^^ cx 
+ 0^ — 2a(3 sm{20). Their contribution to the relax- 
ation rate is enhanced for fields along [110], especially as 
(3 approaches a for a fixed l^o (see Fig. [5]). This is easy to 
understand by looking back at Fig. [D^b) . The relaxation 
of the spin occurs due to dot displacements along direc- 
tions that flip the spin, which for spins along the [110] 
direction (6* = 37r/4) are displacements along [110] itself 
(recall the geometric interpretation as a rolling sphere). 
As we see from Fig. HJb) such angles are the most ef- 
fective to induce spin flips, so the relaxation and also 
decoherence rates increase for fields in those directions, 
especially for the highly anisotropic case a p. In con- 
trast, as is obvious from Eq. p7|) and Fig. [S^b), at low 
fields the rates are dominated by the geometric term that 
is independent of 9. 



V. ELECTRIC SPIN MANIPULATION 

In the previous sections we investigated how in the 
presence of spin-orbit interaction stochastic and fluctu- 
ating fields lead to decay and decoherence of electron spin 
states. In the present section we discuss how one can use 
this effect in a constructive way to control spins purely 
by electric fields. For this purpose one should displace 
the quantum dot which confines the electron. Unfor- 
tunately, under realistic conditions such displacements 
are rather small compared to I so- Nevertheless, a series 
of small closed paths can be designed to take the spin 
of the confined electron to an arbitrary final state. In 
this way one can realize an all-electrical universal single 
qubit gate, under the realistic condition that the scat- 
tering mean free path in the 2DEG (typically in the ^m 
range) is much larger than the confinement lengthscale 
xq. Another option that we shall discuss in this section 
is to move an electron in a system of quantum dots con- 
trolled by gate voltages, and in this way manipulate its 
spin state. 

Closed trajectories of the electron (in a confining dot) 
covering an area A induce a spin precession approxi- 
mately given by Eq. ((29|) around the z axis, in the same 
way as if a constant magnetic field was applied in this di- 
rection. In order to induce arbitrary spin rotations, e.g. 
spin flips, we have to rely on more complicated paths. 
An example is a path composed of sum of closed loops of 
period Tf and area Af and another, much slower, closed 
trajectory of period T, and area As, i.e., a spirograph- 
type of path. One can demonstrate that by properly 
choosing the relation between frequencies and trajecto- 
ries, after a long enough driving time the spin can be 
driven to an arbitrary final state. The optimal relation 
for a spin-flip operation is Tf/Ts = 2Af / {irl'^^). Due 



to the adiabaticity requirement, however, the minimal 
spin-flip time using this method for realistic values of the 
maximum displacement becomes several orders of mag- 
nitude slower than current flip times achieved using ESR 
techniques'^ . The operation time can be reduced by using 
hctcrostructures with larger spin-orbit couplings, such as 
InAs. It can also be reduced substantially if we substitute 
the effect of the fast path component by an equivalent ex- 
ternal static magnetic field along the z direction, in which 
case the technique resembles closely previously proposed 
ac electric-field generalizations of ESR techniques^. 

One can, however, also use different methods to trans- 
port a confined electron over distances comparable or 
greater than the spin-orbit length Igo ~ 3/im. Surface 
acoustic waves, e.g., have been used to move electrons 
over large distances and to rotate their spins^'*. Another 
possibility involves a ring of several tunnel-coupled quan- 
tum dots spanning distances comparable to Iso- The elec- 
tron can then be adiabatically shifted around the ring by 
appropriate time-dependent gate voltages. As we show 
below, such a manipulation can result in a large spin ro- 
tation, and thus provide a completely different principle 
than ESR, since no resonant AC fields would be involved. 
For this reason it could be expected to enable relatively 
fast coherent spin manipulation, since, unlike ESR, a sin- 
gle pumping cycle could be enough to induce a complete 
spin flip. 

The discussion of pumping around the ring of dots can 
be reduced to sequential pumping processes between ad- 
jacent dots. We shall therefore analyze the precession 
of the (pseudo)spin of an electron that is transferred be- 
tween two dots of sizes and separation smaller than ^so- 
We further assume that a strong barrier between the dots 
remains present at all times during the pumping process, 
and use a tight-binding approximation. In this spirit, we 
write the Hamiltonian of the double dot structure as 

^W = |^+ E ^4r) + ^i3A,-i<5- + Uext(r,i) , (38) 

a=L,R 

where Va{r) = V{f — r^) denotes the confining poten- 
tial of the two dots at positions Tq and "Tq, and the 
potential Vcxt is generated by the external gate voltages, 
assumed to be small. Within the tight binding approx- 
imation the low energy states of the double dot system 
are a linear combination of the ground state wavefunc- 
tions IV'oi)} of the isolated dots described by the 

Hamiltonians 

H,;n^^ + V,^, + ^pX-^a (39) 

In the considered limit the SO coupling can be treated 
perturbatively, and the Hamiltonians can be par- 

tially diagonalized by the unitary transformation Zj^ /ji — 
(AMl/r^ with 

= -(r-r^)A-i<T 

Z+Ho^Z^ = l^ + V{r-r",)--^L^S. + 0\X-;:r). 
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Here Iso = \/| det Aso| and is the angular momentum 
with respect to the center of the dot. Let us now assume 
that the confining potentials Va are cylindrically sym- 
metrical and also that the ground states of Ha are in the 
Lz = sector. Then the ground state of the individual 
dots can be approximated, to order as 



(40) 



in terms of the spinor jcr) and the orbital eigenstate |^o) 
of Ha with the spin-orbit coupling set to zero. This al- 
lows us to evaluate the 4x4 matrix for the truncated 
H{t) and express it as 



Hit) = 



Hll Hlb 
Hrl Hjiii 



(41) 



with the submatrices given by {Hafi)(ja' — 
(0Q^|_ff(i)|0Q^,). Apart from a trivial overall shift 
of the energy, the diagonal blocks can be written as 



H 



LL 



'H 



BR 



2 



1 
1 



(42) 



v{t) being the potential difference between the dots. 
However, spin-orbit coupling generates a spin texture 
for the dot eigenstates, and results in a non-trivial spin- 
mixing in the in the hopping submatrix, a ^ a': 



H. 



LR 



= H+^ = {^t'\mM. 



{^^,\H{m^){a'\ 



; Sf x~ 



(43) 
(44) 



where 5r = tr — tl is the vector connecting the two 
dots. In the second line we used the expression (|40)) and 
exploited the fact that the integrals pick up their major 
contributions from the regions r sa and r tr. We 
thus obtain, 



H 



LR 



A{t) ( cos(p) ie-^'^sin(p) 
ie"^ sin(/9) cos(p) 



where p = |A~^ ■ 5r\ is essentially the tunneling dis- 
tance in units of the spin-orbit length and A(t) the spin- 
independent hopping integral (^q | . . . Ivf^) above. The 
angle (j) characterizes the hopping direction and is de- 
fined through the relation 5r ■ — p(cos cj), sin cj)). The 

SU(2) operator e~"''"'^=° oc Hrl is the same we ob- 
tained in Eq. (|22p for the geometric spin precession along 
a straight path connecting the two dots. 

Diagonalizing this Hamiltonian matrix we obtain the 
instantaneous eigenstates {|0„±(i))} with n = and n = 
1 corresponding to the ground and first excited doublets 
of the double dot structure. They are pairwise degenerate 
at any time and have energies eo.cr — —Se{t)/2 and ei,o- — 
Se{t)/2, with Se{t) the splitting between the two doublets, 



5e{t) = ^A{tr+vity 



(45) 



We shall now study the electron spin's evolution 
within the adiabatic approximation: we look for a so- 
lution of the Schrodinger equation in the form \4>{t)) = 



J2n (7 '^■n',rT{t)\(t>na-(t)) . Then the wave function ampli- 
tudes satisfy the equation of motion 

(46) 

7i' ,a' 



I ar. 



+ i A 



+ i V 



dv 



-\4'n'a') 



dA 



If the time-derivatives in this expression are small com- 
pared to the splitting Se the evolution of the confined 
electron is adiabatic and is confined to the lowest dou- 
blet of the double dot. Remarkably, it is possible to write 
down such instantaneous ground states, satisfying the 
conditions (^\6n„,) = (^1^.^) = 0, 



9A 



MA,v) 



(t>o-{A,v) 



~ \ dv 
1 



v/A2 + {Se - vy 

X ((5e — V, 0, — A cos p, Aie*"^ sin p) 
1 



(47) 



(48) 



^A^ + {Se-v)^ 
X (O, (5e — u, Aie"'*^ sinp, — Acosp) . 

With this choice of basis the time evolution is trivial in 
the adiabatic approximation, and apart from an overall 
phase the wave function is simply given by 

m)) ^ a^\^o-iA{t),vm +a+\(j)o+iAit),vm . 

(49) 

Now imagine making an adiabatic sweep, with the po- 
tential difference v going from v = — oo at time t — —oo 
to w = oo at time t — oo. The above states have 
been chosen so that they satisfy the initial condition 
\<Po±{~oo)) — |0o±) at i = — oo, and describe an elec- 
tron localized in the left potential well with (pseudo)spins 
tJ = ±. According to the above expressions, at time 
t = oo, i.e. after the adiabatic potential sweep v oo, 
the electron will be found fully localized in the right dot 
and in the following spin superposition 



/>^j0(-oo)) (50) 



independently of A{t). In other words, the spin un- 
dergoes a spin precession identical to that obtained for 
-5 = upon adiabatically displacing a parabolic con- 
fining potential a distance Sr along a straight line, Eq. 
((22)) . Although in our discussion we assumed a cylindri- 
cal symmetry for the dots, a slightly modified version of 
this discussion carries over to the case of non-cylindrical 
potentials with a spin rotation of comparable size. In 
general, however, the spin rotation is not given by the 
simple expression ([50)) . 

The most straightforward way to build a tunable gate 
on this concept would be to add a magnetic flux across 
a ring of quantum dots as a tuning parameter. Alterna- 
tively one could use backgates on top of the heterostruc- 
ture to tune the Rashba spin-orbit coupling strengths 
around the ring. 
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VI. CONCLUSION 

In this article we have shown that for single electrons 
confined in quantum dots in a 2DEG, which is shifted 
adiabatically along a path by applied or fluctuating elec- 
tric fields, the spin-orbit interaction induces pseudospin 
precession within the ground state Kramer's doublet. In 
the absence of external magnetic fields, the precession de- 
pends solely on the geometrical shape of the trajectory 
of the confined electron. This accumulated non-Abelian 
phase has marked consequences for the spin relaxation 
and decoherence due to electric field fluctuations. In par- 
ticular it leads to a saturation of the relaxation rates 
at vanishing magnetic fields. We have analyzed how 
the properties and power spectrum of the electromag- 
netic fluctuations influences the spin relaxation rates. 
We characterized two different spin decay regimes, domi- 
nated by Ohmic or phonon-induced electric fluctuations, 
respectively, and their crossovers as a function of external 
magnetic field and temperature. 



The geometric spin precession analyzed in this paper 
can also be used to manipulate the spin state purely 
by controlling electric fields. We have shown that arbi- 
trary rotations can be achieved by moving the dot along 
suitable trajectories. To speed-up the process we sug- 
gest moving the electron by adiabatic tunneling between 
quantum dots in multi-dot devices. 
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APPENDIX A: VARIOUS MATRIX ELEMENT 
RELATIONS FOR A PARABOLIC POTENTIAL 

For the parabolic quantum dot Hamiltonian, Eq. (|13p . 
we have the following identity: 



where r is the electron position operator, and Rc{t)) is the dot displacement, assumed to be zero at i = 0. Here Enr 
is the eigenvalue of the instantaneous eigenstates, |T„(t)), which are related to the eigenstates of 7i(0) as |T„(t)) = 
W{t)\Tn). Similarly, the expectation value of r — Rc{t) reads 



{E^-r' - E,,r) (T„W|f - Rcit)K,it)) = (t„(0| [r - Rcit),Hit)\ |t;,(0) = -(t„WIp + A.-^ • aK{t)) 
Combining both equations above we obtain 



1 - (Enr - S„'r')V^O 



A„„"^ • (tJo-It', 



2 /, ,2 '*so \ ' n I 



(Al) 



We made use of the fact that W(t) commutes both with p and &, and therefore (T„(i)|d'|r^, (t)) = (t„|<t|t^,) and 
(r„(t)|p|<,(i)) = (r„|p|<,). ^ 



Specifically, for the ground state doublet (or lowest lying 
two states) this equation reduces to the useful identity 



PappQ = xj 



(A2) 



where t^-t' are usual Pauli matrices. All spin-dressing 
corrections are contained in the matrix A„„ , whose def- 



inition trivially follows from Eq. (|A2p . and simplifies to 
Eq. (pij) in the absence of magnetic field. 



electromagnetic fields are produced by a quantum bath, 
which we shall treat within the path integral formalism. 
We assume that the bath is governed by Hamiltonian 
Ti-B, while the coupling between the bath and quantum 
dot is of dipolar form, V — eE ■ r. We further allow for 
a Zeeman field B coupled to the dot, 



Hi 



P 

2m 



(Bl) 



APPENDIX B: PSEUDOSPIN EVOLUTION 
UNDER QUANTUM FIELDS 



To combine the path integral approach with the adia- 
batic approximation, we write the evolution operator of 
the coupled system as 



In this appendix we generalize Eq. obtained for 

an electron confined in a parabolic well to the case where 



U 



HB+V)t ^ g--i(WD+WB+V)At 



t/At 



(B2) 
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and insert the identity operator 

i = ^VV|£;,a)iD(-E,a|VV^ 



(B3) 



E,a 



after each time shoe. Here Id denotes the identity oper- 
ator acting on the dot, and VV = exp(— ip • Rc)- Note 
that the VV = exp(— ip • Rq) acts both on the bath and 
on the quantum dot, and now Rq — ~eE/mujQ is an op- 
erator instead of a c-number. In Eq. (|B3|) we also made 
use of the fact that the electric field operator E is Her- 
mitian and one can construct a complete basis from its 
eigenstates. However, E being a local operator, every 
such state is infinitely degenerate. We keep track of this 
internal degeneracy by the label a. 

The matrix elements connecting two consecutive iden- 
tity insertions labeled by n -I- 1 and n, take the following 

(£;„+!«„+! I VV;+e-'(^-+^-+^)^*VV;|£;„a„) (B4) 
= (£;„+ia„+i|e-^^«^*|£;„a„) x 



where the operator We^ = e~'P , Rq — —eEn/mujQ 
acts now only on the dots subspace, just as in Sec. lHI ATI 
Using En+i ~ En + AtEn, the second term can be 
expanded in At and written as 



i{HT,+eE„r)Aty^^ 

Hd-P-Rc^ 



e'El 



At 



After re-exponentiating this expression Eq. (|B4p simply 
becomes 



(B5) 



where H'^ = Hs — e^ E^/ {2mLj'^) is effective bath Hamil- 
tonian. 

The rest of the derivation follows the standard con- 
struction of the path integral except that we also insert 
the identity operator (jB3[) before and after the evolution 
operator U (t) . The evolution operator of the dot for fixed 
initial and final bath states Ei and E / , takes the form 



{Ef\U{t)\E,) = / I?[£;]e-^^« (B6) 

J Ei 



numbers a{t) is already integrated out. Finally, the op- 
erator r{e~' ' } inside the integral is acting only on the 
quantum dot, and has exactly the same form as Eq. (fTS]) 
for a classical field. 

Since electric field fluctuations are assumed to be slow 
and the magnetic field is small, the time ordered operator 
in Eq. (|B6p can be approximated by its adiabatic form. 
Projecting this operator to the ground state subspace of 
H with S = 0, we have 

Hv-P- Rc{t') ^ PonvPo - PopPo ■ Rc{t'). (B7) 



After using identity (|A2p of the previous appendix this 
expression reduces to Eq. pop . given in the main text. 



APPENDIX C: DIAGRAMMATIC 
CALCULATION OF RELAXATION AND 
DECOHERENCE RATES 

We can generalize the calculation of the previous ap- 
pendix to compute the evolution of the center of mass 
reduced density matrix, defined as in Eq. p3p . 

PD{t)rr' = (rolTrs [>V+p(t)>v] K) . (CI) 

We use the Schrodinger representation and insert the 
identity operator (|B3p for both the forward and back- 
ward propagation in the expression of the full density 
matrix p{t) = e~'^ *p(0)e*^ *. This leads to the follow- 
ing expression [recall definitions (j3ip and 



where p{Q) = >V~'"p(0)>V denotes the initial center of 
mass density matrix. The integration should be per- 
formed on the complex contour, z = {0 — > t — > 0}, and 
Tpc denotes the time ordering along this contour. 

We define an operator describing the propagation of 
Pd between different times 



Poit)-, 



E 



n{t.^)r,r,^r^rCPD{Q)r^r-. (C2) 



In perturbation theory in the geometric coupling Hq be- 
tween the dot and the electric field we can construct a 
corresponding Dyson equation for the propagator H, 



We,t {e-/o''*'[«-p«c(0]}w+ . H(i,o) = n°(t,o)-f / dt,dt,n\tM)ntiMmh,o) 

Jo 



Here the functional integral is performed over all possible 
paths of bath states with definite E{t) compatible with 
the endpoints Ei and Ef, each of them corresponding 
to a different path C of the displacement Rc{t). Bath 
paths begin a,t t' — and end at t' = t. The weight 
e-'^B = e-'SB[^?(*)] = /I?[a(t)]e-*^B comes 
from the prefactor in (jBSp , and involves the effective bath 
Hamiltonian H'g. Its dependence on the extra quantum 



where 



H"(t,0) 



-itHz 



is the bare propagator. In the following, we shall assume 
that the initial center of mass density matrix factorizes as 
p(0) = /5d(0) X pbath{0), where pbath{0) represents some 
density matrix of a non-interacting Gaussian heat bath. 
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FIG. 7: Diagrams involved in the evaluation of the relaxation 
times. The 4x4 self-energy matrix E that dresses the propa- 
gator n is calculated to fourth order in the dot-bath coupling 
(blue dots). Contractions of bath fields in T,^*'> can be classi- 
fied in two distinct types (upper and lower rows of diagrams). 
Times {0, T\,T2,t} are ordered, and n, t2 must be integrated. 



Then the 'self-energy' T,{t) has an expansion in Feyn- 
man diagrams along the K contour, with vertices at 
branch s = ± corresponding to terms —s iHg, in the ex- 
pansion (see Fig. [7]). Differentiating with respect to t and 
defining the Liouvillian as Lopoit) = -^Ilo{t)p]j{0) = 
i [pnit), Hz], one arrives at the master equation 



poit) = Lopoit) + j:{t~t')pD{t') . 
Jo 

This equation can be simplified under the Markovian 
approximation^, where we assume that relaxation and 
dephasing are slow and therefore replace poit') ~ 
gio(t ^*-)pj-i{t) in the second integrand to yield 

Poit) = Lopoit) + Tpoit) , 
with the Bloch-Redfield tensor defined as 

/•OO 

r / j:{t)e-^°* dt . 

Jo 

Relaxation and dephasing times Ti and T2 are then triv- 
ially related to this tensor F as^ 



t; 



^2 - -R-c ^ti^n 



(C3) 
(C4) 



A 'rotating wave approximation' is implicit in these rela- 
tions that requires T^^, T^^ < ujb in our particular case, 
see Sec. HTCl 

Let us now perform a calculation of Ti and T2 to fourth 
order in the coupling between the dot and the electromag- 
netic field, for a small quantum dot with negligible spin 
dressing (Ago — \o and B — B). Assuming an in-plane 
magnetic field at an angle with respect to the direction 
[100] we have 



Hz = ^M-B-B (cos 6* CTa; + sine dy) 



The eigenvalues of Hz are ±wb/2, where ujb = —giJ-sB. 
In the particular doublet basis which diagonalizes Hz we 
can write Hz = —TzWb/^, and Xl" equals 



with the indices ordered as {ttj Tii ih il})- This basis is 
rotated with respect to the original one by an operator 
Un — exp [— ^(sin0 ~ cos 9 ^y)], so that our vertex 
—isHq in branch s = ±, reads in this basis 




-isHc = -is ^ Rcf,{ts)Ko^^MR^i'U^ 

^.v—x,y 



tJ.=x,y 



xo 



where 7^ = {7;^, 7^ , 7^} = {- sin6i7^, cos 6*7;;-, 7ji} and 
the relevant transverse and parallel dimensionless cou- 
plings, 7^ and 7^, are defined by 

— a cos 6' + /3 sin 6* /3 cos — a sin 6* 
— /3cos0 — asinS a cos + a sin 

(C5) 

For compactness, we shall use the vector notation 7^'ll = 
Klx''^ tly'^^ in the rest of this appendix. The corre- 
sponding vertex matrices V+(i) = —iHQ® \ and V_(t) = 
i\® Hq at time t at branch s — ± that enter the ex- 
pansion of Yj{t) (are denoted by blue dots in Fig. [7l and 
read 



/ 




e^S;^ 





\ 







-ill 


e-'lil 








111 







V 





-e-'S^ 


^ll 




/ 


ill e 







\ 




-e^'ll^ 


-111 












i^l e 


-7f 




V 





-e*''7^ 





In this notation wc have the following relation for the 
n-th order contribution to Bloch-Redfield tensor F in the 
dot-bath coupling V{t) = V+{t) + V-{t) 



p(") 



. dtr, 



(C6) 



(y(t„)nO(i„ - t„_i) . . . y(o)n°(o - i„)>j 



Here tn > ■ ■ • > ^2, i-e. a time-ordered integral over the 
internal times is implicit, and the bracket {■ ■ ■)b denotes 
averaging over the bath fields Rc^n i-e. pairwise con- 
tractions for a non-interacting (Gaussian) bath. Only 
connected contributions must be taken into account in 
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K,(G),,a)2) 



Kj(o),,(a2] 









-1 1 
(D,/C0e 



FIG. 8: The kernels K\ and K2 involved in the result for T^'^'j 1 
Eq. (|C10|) for generic values of the coupling constants 7^'". 
A resolving imaginary part e = 0.0 la's was used to make the 
delta functions visible. 



the averaging, and all odd n contributions average to 
zero. The n = 2 and n = A diagrams are represented 
in Fig. [71 The n — 4 case has two types of distinct 
contractions that correspond to contracting the vertices 
at ti = and with each-other and with internal ver- 
tices, respectively (represented in two separate rows in 
Fig. [7]). Each contraction gives a non- interacting bath 
Green's function defined in the Heisenberg picture as 



TKRcf.it'^.)RcAis)PB{Q) 
V[E]e-^^'^Rc^it'^,)Rc.its) 



Here once more s, s' = ± refer to the time branch, and 

the time evolution of the operators Rcfi{t) in the first 
line is governed by H'g. For simplicity, we assume an 
isotropic bath so that G^'^ — ^^^G"''*. 

In the literature, one commonly denotes in the off- 
diagonal components as G"^ = G"' and G^ = G ^. 
The spectral function is related to these as A = i(G^ — 
G^), while the Keldysh propagator is given by = 
G^ -I- G< . Since E{t) is a real bosonic field, in equilib- 
rium we have 



G>H = G<i-uj)^-iil + nBiLo))A^iLo), 



G^{uj) = -icoth 



2kBT 



(C7) 



where ns(aj) is the Bose distribution, and Ajfiu) is the 
spectral function of the rescaled bath field R/xq. This 
spectral function can be easily related to the spectral 
function A{u!) of the dimensionless electric field, exoE/ojQ 
as A^{lj) = lj^A{lj). Both spectral functions are odd in 

UJ. 

The second order results for Ti and T2 take the known 



form^^ 



(2)- 



2|7^|^*G^(^b) 



Tr 



(C8) 
(C9) 



The fourth order contributions to F are obtained by sum- 
ming all corresponding diagrams. They are rather in- 
volved, but have the general structure 



(4)- 

1:2 



(CIO) 

The kernel contains delta functions of wi , 0-12 . In Fig. 
[8] we plot the lines along which these delta functions pick 
up their contribution for relaxation and dephasing. We 
can distinguish two types of lines. The diagonal lines in 
the 2nd and 4th quadrants lead to convergent integrals, 
since G-^ {uj)G^ {—lu) goes exponentially to zero at large 
values of to. Diagonal lines with positive slope and also 
horizontal and vertical lines across the origin cancel out 
exactly in Ki. The remaining horizontal and vertical 
lines exhibit ultraviolet divergencies, and give therefore 
cutoff dependent prefactors that multiply G^{ztuj) and 
G^(0). These terms can thus be reabsorbed in the second 
order result by simply renormalizing the couplings •y-^'K 
The final result can then be expressed in terms of these 
renormalized constants as 



(Cll) 



-^2 



-217^1^17" l^iFf (lob) + 2 (7^ • 7" j ^Fi' (lob) 
:Ti-i + |7"|'^G^(0) + h^l^tFiuB) (C12) 



where the F functions are given by the following convo- 
lutions of G^(cj) 

J^^2M = F>2(C.) + F>2(-W), 

Re 







1 






(f 


1- - 


f io) 



.>M^ / |^.c> (I (I 



Re 



1 



1 



^-uj' -iO^+uj' + iO 



F{u;) 



"2^ 



iG> {-uj')iG> {uj') 



Re 



1 



1 



CO — Lv' — iO u! + oj' + iO 



.(C13) 



The Green's functions G^ and G^ can be expressed in 
terms of the spectral function A{uj) of the dimensionless 
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electric field exoE/ujo, introduced below Eq. (|C7 



It can be checked that and Ey are indeed independent 
in this model. 



iG'^Lu) = coth 



Using (^(t)6+(0)) 



\1 



2kBT 

Inserting these last equations into the expression (|C13|) 
one arrives at the results in Eqs. ([55)1 and The 
spectral function A(uj) for electromagnetic fluctuations 
generated by piezoelectric phonons is computed in the 
next appendix. 

APPENDIX D: PHONON BATH PROPERTIES 

Let us consider the fluctuating electric field induced 
by the phonons in the sample holding the quantum dot. 
The electric field acting on the confined electron is the 
gradient of the potential generated by these phonons, 
eEfj, = —\/fj,UP^{x,y). It has two contributions, one 
from a longitudinal mode and another from two trans- 
verse modes 



uj^A{u;) . (C14) {b+{t)bq{0)) = e'"''*ni3(wg) and = viq we see 



(Dl) 



with A = l,ti,t2 indicating the mode. The coupling to 
the longitudinal is given by^ 



o \ 3 I Jil^dw) 
2pviq \ q-' J 

while the coupling to the transverse mode is 



iqxQyqz 



'14 



Jiqzdu 



2pvtq 



Here q = \q\, V is volume, is the depth of the quan- 
tum well where the 2-dimensional electron gas is con- 
fined. The function J{x) — 0(1 — x'^), 8 being the 
Heaviside function, qualitatively accounts for the trun- 
cation of the phonon spectrum out of the 2DEG plane^^i 
The physical origin of this cut-off is that phonons hav- 
ing a wave- vector component larger than ~ l/du, along 
the z direction cannot couple efficiently to the confined 
electron, since their wave function oscillates too quickly. 
The density g, the sound velocities vi and Vt, and the 
piezoelectric constant hi4 in the expressions above are 
material-dependent parameters, which depend on the 
particular heterostructure used to define the quantum 
dot: for a typical GaAs/AlGaAs heterostructures p — 
5.3 • 10^ Kg/m3, VI = 4.73 • 10^ m/s, = 3.35 • 10^ m/s, 
and hi4 = 1.4 • 10^ V/m. 

The noise power S>{t'~t) = (E^{t'_)E^{t+))Be'^xl/ujl 
of the normalized electric field can be expressed as 

S>{t) 



™2 
-''0 



(&+AW&q,A(0)) 



(N,aW&+a(0)> 



that the only dependence on the orientation of the wave 
vector vector q appears through M. We can therefore 
introduce spherical coordinates in the continuum limit 
and integrate with respect to the angular variables. For 
the longitudinal phonons this yields 



ujgpvi 2107r 



where the function ji{qdw) above comes from the trun- 
cation of the phonon spectrum, and is given by 



]i{x > 1) = 



-35 4- 135a;2 - 189x'' + lO^x^ 



16x9 



3i{0 < X < 1) = 1 



(D3) 



The two transverse modes give a similar contribution 



sf{t) 



ujQpVt 210 



-r^Jt{qd^)q' (D4) 



X (e^"-**nBKt) + e-^--'*[l + 71bM]) , 



with 



jt{x > 1) = 



105 - 300x2 _^ 294^4 - 140x6 105x^ 



64x9 



jt(0 < X < 1) = 1 



(D5) 



The cut-off functions ji (x) and jt (x) play a role similar to 



0(1 — X ), but have algebraic tails: ji{x) 



for large 



X, while jt{x) ~ 1/x. The phonon spectrum is defined as 
usual by uiq^x = v\q. 

Taking the Fourier transform of the correlation func- 
tions above, we can identify the spectral function Api^^Lu), 



S>{ij) = [l + nB{ij)]A,M 
Aph{uj) 



ph 



2107rp 



(D6) 
(D7) 



01 



where 10^,1 = vi/dw and uj^^t = Vt/dw For frequen- 
cies smaller than the cutoffs lo^^i we obtain Aph = 2.5 • 
10"^K~2nm~2 for GaAs/AlGaAs quantum wells. Note 
that in this calculation we neglected the lateral size of the 
dot. That provides an additional cut-off for the in-plane 
phonon momentum for larger values of the frequency lo. 
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